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Abstract 



The Glauber dynamics of various models (REM-like trap models, Brownian motion, 
BM model, Ising chain and SK model) is analyzed in relation with the existence of ageing. 
From a finite size Glauber matrix, we calculate a time t w (N) after which the system has 
relaxed to the equilibrium state. The case of metastability is also discussed. If the only 
■ non zero overlaps between pure states are only self-overlaps (REM-like trap models, BM 

model), the existence or absence of ageing depends only on the behavior of the density of 
eigenvalues for small eigenvalues. We have carried out a detailed numerical and analytical 
analysis of the density of eigenvalues of the REM-like trap models. In this case, we show 
that the behavior of the density of eigenvalues for typical trap realizations is related to 
the spectral dimension of the equivalent random walk model. 
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1 Introduction 



Various approaches have been used to understand ageing experiments || ||] @: droplet picture §, 
mean field model § f§, and trap models @ |] § @. A geing is by definition the property that 
the two points correlation function 

C(t w ,t) = {0(t w + t)0(t w )) (1) 

depends explicitly on the waiting time t w . O is a macroscopic observable, for instance the magneti- 
zation in the zero field cooled experiment |Q or the thermoremanent magnetization experiment Q. 
On the theoretical level, dynamics can be implemented in several ways. One possibility is Langevin 
dynamics which was used for instance in [11] where mean field ageing was first obtained, or in [12] 



where ageing in the presence of flat directions was studied. Another possibility is Glauber dynamics 
[ |13| . Within the framework of Glauber dynamics and given the probability distribution p({ej},0), 
for the system to be in any microcanonical state at time t = 0, the probability distribution at time 
t is given by 

p({e i },t)=exp(Tt)p({e i },0), (2) 

where T is the Glauber matrix. By "microcanonical states", we mean a state such as the system is 
in a given microcanonical configuration with probability one. 

This article is devoted to the study of the Glauber dynamics of various models, phenomenological 
or microscopic, in connection with ageing. More precisely, for a finite size system of size N, we 
calculate from the Glauber matrix a time t w (N) such as the finite size system equilibrates on a time 
of the order of t w (N), starting from out-of-equilibrium initial conditions at time t = 0. In terms of 
ageing, the correlation function (|l|) depends explicitely on the waiting time t w provided the system 
is out-of-equilibrium at time t w , namely t w is smaller than t w (N). In other words, ageing for a 
finite size system is interupted after t w (N) [see JtJ where the denomination "interrupted ageing" 
was proposed]. In order to understand ageing in the thermodynamic limit (N — > oo), the variations 
of t w (N) with the system size have to be carefully studied. If t w (N) goes to a constant in the 
thermodynamic limit, ageing is interrupted, whereas if t w (N) diverges ageing is not interrupted in 
the thermodynamic limit: in the thermodynamic limit, the system reaches equilibrium on infinite 
time scales. We also discuss the possibility of metastability without ageing on experimental time 
scales, as it is the case for instance in diamond. 

In general, t w (N) depends on the eigenvalues of the Glauber matrix, but also on the overlaps 
between microcanonical states, and the transformation matrix elements between the eigenstates of 
the Glauber matrix and the microcanonical states. However, in some simple cases, with only self- 
overlaping states, the time t w (N) can be expressed only in terms of the density of eigenvalues of the 
Glauber matrix. 
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The form of the Glauber matrix will be given in each case in the core of the paper. Whatever the 
system, the eigenvalues of T are negative. The spectrum is bounded above by zero (zero is always 
an eigenvalue, corresponding to the equilibrium Boltzmann distribution). The eigenvalues of the 
Glauber matrix have the dimension of the inverse of a time, so that they can be viewed as the inverse 
relaxation times of the system. Since we are interested in long time behavior, a special attention will 
be paid to the smallest eigenvalues. 

This article is organized as follows. In section & we present a detailed study of the average density 
of small eigenvalues of the random energy model (REM) like trap model proposed and studied in 
0, |, | and inspired by the REM model 0. 

In the case of the REM-like trap model, two regimes are predicted and checked numerically for 
the density of small eigenvalues in the glass phase: an Arrhenius-like regime and a non Arrhenius 
regime for the smallest eigenvalues. By Arrhenius-like regime, we mean a regime in which the density 
of eigenvalues is the same as the density of eigenvalues calculated from the Arrhenius law, assuming 
a correspondence between the eigenvalues of the Glauber matrix and the inverse relaxation times 



of the system. The eigenvalue density in the Arrhenius regime are given in section 2.1. We also 
establish a connection between the behavior of the density of small eigenvalues of REM-like trap 
models and the spectral dimension of the associated diffusion process. 

Next, in section ||[ we give a general relation for the existence of ageing from a Glauber dynamics 
point of view. This section is the main part of the present article. When applied to the REM-like 
trap model, our criterium is equivalent to the fact noticed in Q that the average trapping time 
diverges below T g . The rest of the paper is devoted to examine the case of other models: Brownian 
motion which was shown in []Fj| to be an ageing phenomenon, a model introduced by Barrat and 
Mezard (BM model), the one dimensional Ising chain (in this case, our results are consistent 
with the existence of interrupted ageing), and the Sherrington-Kirkpatrick model (SK model). In all 
these cases, we find consistent results between the existence of ageing calculated from the behavior of 
the two times autocorrelation functions and our criterium based on the Glauber matrix. We end-up 
with some final remarks. 



2 Density of small eigenvalues of trap models (REM-like trap mod- 
els and trap models with other topologies) 

Trap models can be defined without a priori reference to a microscopic model. However, some trap 
models, such as the REM-like trap model § are related to the Parisi solution of spin glasses: even 
though these models are toy models, some of them retain some feature of the dynamics of microscopic 
models. The REM-like trap model consists of a trap model in infinite dimension, that is, the system 
can "hop" from one trap to any other trap, according to the rules of the dynamics. Other toplogies 
will be analyzed below, in reduced dimension. The REM-like trap model corresponds to a single 
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replica symmetry breaking step. Generalizations of the REM-like trap model including multi replica 
symmetry breaking steps were studied in || and their relevance to explain experimental data was 
shown. 

The models of this section consist of N traps of depth Ei < 0, with the Ei chosen among the 
distribution P(E) = exp(E). The reason why this particular energy trap distribution is chosen is 
that there exists a finite transition temperature between a high and a low temperature glassy phase 
[ 10 1 . The dynamics rules are given by Tjj oc exp {(3Ej) for i ^ j if the traps i and j are connected 
and zero otherwise. The diagonal terms are defined by 

N 

E T M = °' (3) 

k=l 

which enforces probability conservation. Equation (|3|) is a general constraint for any Glauber dy- 
namics. Another constraint is the detailed balance condition 

T hi exp {-(3Ei) = T id exp (-(3Ej). (4) 

It is clear that the aforementioned rules for the dynamics of the REM-like trap model satisfy detailed 
balance. However, as we shall see below, other types of dynamics also satisfy detailed balance and 
probability conservation. As a consequence of detailed balance, the eigenvectors of the Glauber 
matrix are not physical states (in the sense that they are not probability distributions) except for the 
Boltzmann distribution (the stationary eigenvector corresponding to the A = eigenvalue). Notice 
that if T is a Glauber matrix xT, with x a positive number, is also a Glauber matrix. Multiplying 
T by x amounts to changing the unit time. The proportionality constant is imposed by physical 
considerations. For the models considered in this section, this constant is chosen such as the diagonal 
elements of the Glauber matrix do not scale with the number N of traps in the large N limit. For 
instance, in the infinite dimensional case, we choose Tjj = exp ((3Ej)/N. This dynamics will be 
studied with different phase space topologies: infinite dimension, where all the traps are connected 
(section |2.2| ), one dimensional topology, where a given trap is connected only to its two nearest 
neighbors (section |2.3| ), and also one dimensional topology with "ordered" traps, in the sense that 
the trap energies have been ordered E\ < ... < En, and the trap at site i has an energy E^. 

A useful property that we will use throughout this paper is that Glauber matrices can be 
mapped onto symmetric matrices with the same spectrum. Let us call V the diagonal matrix 
Vij = exp (PEi/2)Sij. Then, because of detailed balance, the matrix T = VTV -1 is symmetric 
and has the same spectrum as T. Using a quantum mechanical analogy, we denote by |ej) the 
microcanonical state corresponding the system in trap i with a probability unity. The analogy to 
quantum mechanics is useful, even though there is no hermitian structure. In particular, any linear 
combination of microcanonical states is not a probability distribution. 
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2.1 Arrhenius regime 

We propose ourselves to carry out a detailed study of the density of eigenvalues of the REM-like trap 
model and related models in reduced dimension. Since the Arrhenius law is introduced "by hand" in 
the dynamics, it is interesting to question in which regime the density of eigenvalues is given by the 
Arrhenius law, assuming a correspondence r = — 1/A between the eigenvalues of the Glauber matrix 
and relaxation times. We will naturally call Arrhenius regime such a regime. As we will see below, 
the density of eigenvalues is not always in the Arrhenius regime. 

In this Arrhenius regime, and by definition, the off-diagonal terms are neglectable as far as the 
density of eigenvalues is concerned. The diagonal coefficients of the Glauber matrix are proportional 
to exp(/3£'j), so that the distribution of eigenvalues is 

P(ln|T|) cxexp(-j-ln|T|). (5) 

P 

This distribution is the same as the one expected from an Arrhenius law with life times Tj oc 
exp (—(3Ei). If one assumes an Arrhenius picture, the distribution of relaxation times is 

P ir) - — Lj, (6) 
which is equivalent to (||) . As noticed in |J , the average trapping time diverges if j3 > 1 , leading to 



ageing. Distribution of trapping times with no first moment were already studied in [18|, and shown 
to lead to anomalous diffusion, namely to a non trivial spectral dimension, even above the upper 
critical dimension. We will compare our results for the REM-like trap model with the results of [18] 
derived of the context of anomalous diffusion on site-disordered hypercubic lattices. 

The density of eigenvalues is thus described by the Arrhenius law provided the off-diagonal terms 
of the Glauber matrix are a small perturbation, so that the density of eigenvalues of the Glauber 
matrix is given by the diagonal terms only. We now examine several cases with different phase 
space topologies and determine in the different cases on which condition the density of eigenvalues 
is determined by the Arrhenius picture. We also compare our analytical approach to numerical 
diagonalizations of the symmetrized matrix T. 

2.2 Infinite dimensional case: REM-like trap model 

In this section, the system can hop from one trap to any other trap (infinite dimension). The diagonal 
coefficients are 

N - 1 

T i4 = = — exp (J3Ei), (7) 

and the off-diagonal coefficients of the symmetrized Glauber matrix are 

r*,; = ^ exp (8) 
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2.2.1 Determination of the Arrhenius regime 

We treat the off-diagonal coefficients T{ j (i ^ j) as a perturbation. This perturbative treatment is 
consistent since we are looking for a regime in which the off-diagonal coefficients do not modify the 
density of eigenvalues with respect to the purely diagonal case. The second order correction to the 
eigenvalue Tj j is 

A2) = 1 V exp WEj + E k )/2) 

* 2N(N-l)^smh((3(E k -E t )/2y [> 

the first order corrections being identically zero since the perturbation consists only of off-diagonal 
elements. We assume that the density of trap energies in the vicinity of Ei is large enough for the 
discrete summation to be replaced by an integral. If this condition is not satisfied, the second order 

(2) 

corrections 5\ may have large fluctuations because the divergence for E k — > Ei is no more cancelled. 
The condition that the unperturbed eigenvalue density is large in the vicinity of Ei simply reads 
N exp (Ei) ^> 1. We assume in the rest of this section that this condition is fulfilled, and we can 
safely replace the discrete sum by an integral: 

A (2) _ 1 f° exp((3(E i + E)/2) E 

°* ~ 2(N - 1) i-oo sinh(/3(£ - ^)/2) 1 j 

The singularities for E = Ef and E = Ey cancel each other, and we can approximate 5^ as 

4° ^ ^^yS exp ((1 + (11) 

with B = 4(1 + (3/2)/ (3 2 . We have used the following approximation: if A is small, then 

exp (Ax) ^ 2\A 
-A sinh (fix) [a 

and we use A = 2/(3 as a cut-off. Clearly ( |TT| ) is only an approximation, but we are only interested 

(2) 

in the dominant behavior of as a function of Ei. Now, we want to compare the approximation 
(0) of Sf } to the typical level spacing between T^j and Tj+i^+i, where we assume that the traps are 
ordered such as E\ < E% < ... < En. To to so, we first calculate the distribution of 

N—l 

\T\ = ^ r exp(PE), (13) 



and we obtain P(\T\) = A^ 1 ^- 1 , with 



1 / N \ 1//3 

^=4h^-T • (14) 



P\N-1. 

Notice here that P(|T|) diverges as \T\ goes to zero if (3 > 1. This accumulation of small relaxation 
times is related to the existence a glass transition for this model for (3 = 1. Let us now calculate 
the level spacing statistics P\^i(S), which is the density probability that the spacing between the 



5 



absolute value of the eigenvalue T and the absolute value of the next eigenvalue is S. Because of the 
statistical independence of the trap depths, Pm (S) is readily obtained as 

In this expression, k is the number of levels in the interval [\T\ + S, \T\ + S + 5<S] (with > 1). 
The other terms correspond to the condition that there is no eigenvalue in the interval [|T|, \T\ + S]. 
Assuming that \T\ is small compared to unity, and that ./V is large, we get 

P m (S)^^N(\T\ + S)^-\ (16) 

The off-diagonal terms are irrelevant to the density of eigenvalues provided 

V{\T\)= / 1 P\ T \{S)dS (17) 
Jo 

is much smaller than unity. This condition reads Be 2Ei / (3 <C 1, with B defined in (11). Since we are 



looking for a necessary condition valid if j3 > 1, we neglect the temperature-dependent prefactor and 
conclude that the density of states is in the Arrhenius regime only for traps such that 1/N -C e E -C 1. 
Notice that the traps such as N exp B < 1 are deeper than the typical deepest trap (of the order 
of — IniV). The Arrhenius picture is thus valid for a typical trap realization, with a lowest trap 
not deeper than — IniV. However, non typical disorder realizations may be generated. For these 
realizations, the lowest eigenvalues are not described by an Arrhenius picture. The existence or 
absence of the Arrhenius regime even for a very large number of traps is thus sample dependent. 
This non self averaging behavior is reminiscent on the sample dependence of the quantity 

N 

y = Y. w l ( 18 ) 



with Wi the probability to find the system in trap i [17]. 
2.2.2 Numerical tests 

The density of small eigenvalues of the REM-like trap model is pictured on figure |l[ For practical 
purposes, we do not calculate the density P(ln|A|) since this quantity is too small for small eigen- 
values (see (||)). We rather calculate Pe(ln|A|) = WB-P(ln | A|), where Wb is the Boltzmann weight 
associated to the eigenvalue A, which is calculated as the average of the Boltzmann weight operator 

N 

W B = J2\^e~ l3Ei (ei\ (19) 
i=i 

over the eigenstate associated to the eigenvalue A. In the Arrhenius picture, we have 

FB(m|A|)oce Xp (Q-l)ln|A|). (20) 
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Since (5 > 1, -Pb (In |A|) increases exponentially as a function of In |A| when A — > 0~. We see on figure 
H| the existence of the Arrhenius regime, in good agreement with the condition 1/N <C e E <C 1. We 
can now fit the density of eigenvalues in the Arrhenius regime to the form 

P B (ln|A|) ocexp(-aln|A|), (21) 

and see if we recover a = 1 — T, as predicted from the Arrhenius law. The fit to the Arrhenius law 
is shown on figure [| The variations of the exponent a as a function of the inverse temperature are 
plotted on figure ||, and we conclude to a good agreement with the Arrhenius law. 

2.2.3 Relation with the anomalous diffusion on site disordered lattices 

Random walks on a finite d-dimensional lattice with a broad distribution of waiting times have been 



analyzed in the past (see [18] and references therein). As we shall see, we can relate our results 
concerning the existence of a Arrhenius-like regime for the REM-like trap model to the results of 
[Q8|, where the long time behavior of random walk on such lattices was investigated by means of 



the real space renormalization method introduced in IS]. The REM-like trap model falls in the 
universality class described by Machta since the ingredients are the same in the two cases, rather 
than the universality class of random walks in a random field force [^]. In order to make contact 
with the work of Machta, we first calculate the probability P (t) of being at the origin: 



A' 



P^ EM \t) = ^£*>(M«,o) = -^E eA ^ (22) 

i=l a 

where V(ei,t\ei, 0) is the probability to find the system in the trap i at time t, starting from the trap 
i at time 0, and {A Q } is the set of eigenvalues of the Glauber matrix. We have shown that, in the 
infinite dimensional case, the typical small eigenvalue density (relevant to the long time behavior of 
Po(t)) has the Arrhenius form p(A) ~ X 7 " 1 . A simple scaling analyzis then shows that 

P^ REM \t)^r T (23) 

in the long time limit. This result was derived for a graph such as any site is connected with any 
other site, that is we started from the begining from an infinite dimensional situation. In order to 
extrapolate this result to a finite dimension d larger than the upper critical dimension 2 (as shown 
in |18||) we first notice that, within the framework of the REM-like trap model, the microscopic time 
ro scales like ~ a d , with a the lattice spacing, whereas, in the context of random walks, 

the correct scaling is Tq RW ^ ~ a 2 [an illustration of the validity of such a scaling will be given in 
section |j. This implies that Tq RW ^ ~ (TQ REM ^) 2 / d , which gives the correct prescription to go from 
the REM-like trap model in infinite dimension to the finite d dimensional random walk model: one 
should replace the time t in p^ ) by t d l 2 . This leads to 

p (RW) ^ t -Td/2 (24) 
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for d larger than the upper critical dimension 2. This is exactly what was found in [18| using a real 
space renormalization group approach, with a non trivial spectral dimension d s = Td. 

2.2.4 Structure of the eigenvectors of the REM like trap model 

We aim to calculate the structure of the eigenstates of the REM-like trap model. To do so, we 
decompose a given microcanonical state |ej) into the eigenvectors and calculate the absolute 

value of the coefficients of this decomposition. The result is plotted on figure |j. We see that a 
microcanonical state |ej) with an eigenvalue Aj is mainly a linear combination of eigenstates with 
eigenvalues close to A« . This observation can be understood in the framework of perturbation theory 
as follows. The first order correction to the eigenvectors are 

|s « = 2(F^T)g S i„h WE I- £i )/2) |e * , • (25) 

so that, at first order 

h) = m ~ whr) g s^{E k -Em^ kY (26) 

At the first order of perturbation theory, the pure states |ej) is thus mainly a linear combination of 
eigenvectors with an eigenvalue close to the eigenvalue of l^}. 

Starting from a pure state |ej) at time t = 0, we can calculate the evolution of the system: 

hit)) = exp(-^-e^t)\ ei ) (27) 



1 x ^ 1 

+ 2(N-l)^ i smh((3(E k -E l )/2) 



exp j - N " 1 e f3El ^j - exp f"- ^^ 1 e? E H 



where we have used (25) and (p6|). As expected from the Arrhenius law, the relaxation time is of the 
order of Tj = exp (—(3Ei), which diverges in the limit (3 — > +oo. We will compare below the structure 
of the eigenvectors of the REM-like trap model to the structure of eigenvectors of another model 
(BM model). 

2.3 One dimensional case with disordered traps 
2.3.1 Absence of Arrhrenius regime 

Assuming periodic boundary conditions, the diagonal terms of the Glauber matrix are given by 

T iti = -2exp{(3E i ). (28) 

The perturbation induced by the off-diagonal terms does not affect significatively the density of eigen- 
values provided the perturbation Tjj = exp ({3(Ei + Ej)/2) is small compared to the level spacing 
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between Ta and 7i.fl,i+i- Since Ej is uncorrelated from Ei, the matrix element Tij must be averaged 
over Ej, so that Tjj is of the order of exp (J3Ei/2), that is 



^'~V 2 ' (29) 
The distribution of eigenvalues of the unperturbed system is 

from what we deduce the level spacing statistics in the long time (small \T\) approximation 

P m (S)^j(\T\+S)VP-\ (31) 
The off-diagonal terms do not affect the density of eigenvalues provided 

P(S)dS < 1, (32) 

o 

which leads to the condition TV 2 exp E <C 1. The average number of traps generated among a 
sequence of N independent traps chosen among the distribution P(E) = expE and fulfilling the 
aforementioned condition is 

/■- 21nJV i 

N L eEdE =N< m 

which is vanishing in the large N limit. We thus conclude that the Arrhenius regime with this 
topology does not exist, that is off-diagonal terms will always modify significatively the density of 
eigenvalues with respect to the Arrhenius case. As we will see below, the Arrhenius regime is replaced 
by a renormalized Arrhenius regime, where the density of eigenvalues is given by (|2l|), but with a 
lower than 1 — T. 

2.3.2 Numerical tests 

Again, we generate numerically the density of small eigenvalues Pb (In |A|). The fit is shown on figure 
^, and the variations of the exponent in the Arrhenius regime are plotted on figure ^. We conclude 
to the existence of a renormalized Arrhenius regime, in the sense that the density of eigenvalues is 
given by ( pip but with an exponent a smaller than 1 — T, indicating a slowing down of the dynamics 
due to the one dimensional topology compared to the infinite dimensional case: in one dimension, to 
go from one trap to a deeper trap, the system has to jump to intermediate traps. 

2.3.3 Relation with the anomalous diffusion on site disordered lattices 



Section 2.2.3 was devoted to the connection between the spectrum of the Glauber matrix in infinite 



dimension and the long time behavior of random walks on finite d dimensional trapped hypercubic 
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lattices above the upper critical dimension. We now analyze the analogy in one dimension. We have 
shown numerically that the density of eigenvalues behaves like p(X) ~ A 1_T " for small A, where T' is 
different from T. In terms of random walks, this result becomes, with the notations of section |2.2.3| 

pl? w \t)~t- T '/\ (34) 

namely the spectral dimension is d s = T' . We now use the results of where d s was calculated: 



2T 

T' = d s = 2v = , (35) 

l + r v ' 



from what we conclude that the coefficient a plotted on figure || is 

l-T 



a 



(36) 



1 + T 

This analytic form is consistent with numerical calculations of the density of eigenvalues. 
2.4 One dimensional case with ordered traps 

We now assume that iV energy traps were generated among the distribution P{E) = expE, and that 
these energies have been ordered in such a way that E\ < E2 < ... < E^. A site k with a trap of 
energy Ej. is connected to the sites k — 1 (with an energy Ek-i) and k + 1 (with an energy Ek+i, 
except for the sites 1 and N which have no neighbor. Intuitively, such a trap ordering should make 
the dynamics faster with respect to the case analyzed in section |2.3| since the system can decrease its 
energy continuously during the relaxation. We have 

Ti.i+i = — \jTi t iTi + i^ + i. (37) 

Since the density of level spacings Pm(S) ( |3l| ) is large for small values of |T|, we expect that Tj+i^+i 
is close to T^j at least if \Ti i\ is small, so that we can use the following approximation: T^, = Tj j/2. 
The off-diagonal terms will not modify the density of eigenvalues with respect to the Arrhenius case 
provided 

r\T\/2 

/ iWS)dS<l, (38) 
Jo 

that is N exp E <C 1. The number of traps fulfilling this condition in a distribution of N traps is on 
average equal to unity, so that only a small number of eigenvalues contribute to the Arrhenius regime. 
Nevertheless, it is still possible to make a statistics over a large ensemble of iV traps distribution and 
to test numerically if the Arrhenius regime can be observed. 

2.4.1 Numerical tests 

We generated numerically the density of small eigenvalues Pb(1u|A|). As shown on figure we still 
can fit -Pb(A) to the form exp (—a In |A|), and the exponent a is close to 1 — T, as in the Arrhenius 
regime, (see figure |3|). We conclude that ordering the traps in one dimension restores the Arrhenius 
regime for the lowest eigenvalues. Again, the existence of the Arrhenius regime is sample dependent. 
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3 General relation between ageing and the Glauber matrix 



We are now going to give a general relation between ageing and the Glauber matrix. This relation is 
a generalization of a criterium given in j|] . In the context of the REM-like trap model, the average 
trapping time 

(t)=t [ dEexp(E(l-p)) (39) 

diverges if T < T g = 1. On the other hand, ageing exists below T g only, in the limit where the 
microscopic time scale To goes to zero. We are now going to argue that a similar criterium can be 
found for any statistical mechanical model, provided the master equation is linear. This is not always 
the case. For instance, the master equation of the Backgammon model in a mean field dynamics at 
zero temperature is not linear 

Let us now start our general argument. The idea is to start with a system in a given microcanonical 
state |ej) and write the Glauber dynamics equation for the evolution of this pure state. First, we call 
P the transformation matrix 

|e*>=E P *A>' ( 4 °) 

a 

where |ej) is a microcanonical state and {|\I' Q ,)} are the eigenstates of the Glauber matrix. By 
microcanonical state, we mean that if the system is in the state |ej), the system is with probability 
one in the microcanonical state e,. In order to avoid confusion, all the labels related to the eigenvectors 
will be denoted by greek letters whereas the labels for pure states will be denoted by latin letters. 
The evolution of the pure state ( |40| ) is readily obtained as 

| ei (t)) = E^ eAa ^1- ( 41 ) 

In other words, the probability to be in microcanonical state ej at time t starting from the micro- 
canonical state ej is 

V( ej , %0) = £ P ija e A «*p-j. (42) 

a 

We now consider the time-dependent correlation function 

= ^E p ( e i>^'°)^' ( 43 ) 

where N is the number of microcanonical states. We use M instead of N to avoid confusion between 
the number of microstates M and the number of sites N in a spin model. In a trap model, N = N, 
whereas in an Ising spin system, N = 2 N . In (|43"1), we have carried out an average over all the 
initial microcanonical states, that is the system is first quenched from an infinite temperature, and 
the correlation ( |43"| ) is averaged over the realizations of the infinite temperature states. In (|43|), qij 
is the overlap between the pure states ei and e^, which has to be defined within each model. As we 
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will see below in section ||, even within a given dynamics, the choice of the overlap depends on what 
type of observables are considered in the two times correlation functions. We now want to calculate 
from (|43| ) a typical time by integrating (^). However, the Boltzmann eigenstate leads to a trivial 
divergence of this integral. We thus need a prescription to eliminate the equilibrium distribution 
in (f43|). We first discuss this prescription in the absence of metastability and postpone for later 
discussion the case of metastability. We define the typical time t w by 



= E P ^Y P ^ (44) 

hi o^Bol. a 

where the summation over the eigenstates excludes the Boltzmann distribution. The overline denotes 
an average over the disorder realizations, if necessary. The significance of t w {N) is clear: for a finite 
size system with N sites, the out-of-equilibrium dynamics occurs before t w {N) and, for times larger 
than t w (N), the finite size system has equilibriated. In other words, t w (N) is the cross-over time 
before which the finite size system exhibits ageing. Notice that in this argument, the system is a 
finite size system. One should thus carefully study finite size effects. Two cases may occur: first, 
t w (N) tends to a constant in the thermodynamic limit. In this case, ageing is interrupted jjj in the 
thermodynamic limit. As we shall see, this is the case for the Ising chain at a finite temperature. 
On the other hand, if t w (N) diverges in the thermodynamic limit, ageing is not interrupted. In this 
case, the exitence of a divergent relaxation time with respect to times and t implies the existence 
of a divergent relaxation time with respect to any finite waiting time t w and t + t w , and thus the 
existence of ageing in the thermodynamic limit. The relation between the existence/absence, the 
nature of ageing in the thermodynamic limit and the behavior of r w (N) was derived on heuristic 
basis. We could not find a rigourous derivation using a general argument. However, in what follows, 
we investigate the behavior of several models (REM-like trap model, random walks, BM model, one 
dimensional Ising chain and SK model) and we find in all these cases that our criterium based on 
the behavior of t w (N) is correct. As far as non interrupted ageing in the thermodynamic limit is 
concerned, we will see below that this is the case for the BM and SK models. This is also the case 
for the REM-like trap model. We have shown in section |2.2| that for all the typical realizations 
of the trap configurations (typical is defined in section |2.2| ), the density of eigenvalues is given by 
the Arrhenius law, which means that, in the case of the REM-like trap model, our criterium for 
the existence of ageing is equivalent to the criterium of Bouchaud M that we have recalled at the 
beginning of this section, provided the overlap between traps is taken to be qij = qEA^i,j, where qEA 
is the Edwards-Anderson order parameter (within the one replica symmetry breaking step picture). 
Notice that the existence of a non Arrhenius regime generated by non typical trap configurations 
does not modify our conclusion regarding the divergence of t^ em (oo) since the non Arrhenius regime 
only occurs for traps such as E < — In N — > — oo if N — > oo. The correlation C^ M (t) is 

'(t)= f° dEe E expf-te^), (45) 



c 



REM/ 
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which is non integrable if (3 < 1. The typical time 



t ™m {n) = _clea ^ J- = - toA /° i piV (A)dA (46) 

a^Boi. 

only depends on the density of eigenvalues PnW for a system of N traps, and is independent on 
the transition matrix, which motivates the study of the density of eigenvalues carried out in section 
||. In general, in any model with only self-overlapping pure states, the existence of ageing is coded 
only in the behavior of the density of small eigenvalues: there is ageing in the thermodynamic limit 
if Poo (A) tends to a constant if A goes to zero (an exemple will be provided in section ||) or if Poo(A) 
diverges if A goes to zero (this is the case for the REM-like trap model already studied in section [2| 

)• 



Let us now discuss the behavior of the correlation function (43) in the presence of metastability. 
We have in mind the case of diamond for instance, which is metastable but in which case there is 
no ageing. We do not have a Glauber matrix description of diamond, but take this case as a generic 
situation for metastability and draw qualitative conclusions. In this case we expect that for a finite 
size system, (|43| ) will first decay to a finite value over a first time scale t£\N). After t$~\n), the 
system reaches its metastable configuration (for instance diamond) and thermalizes in a portion of 
its phase space. We thus expect that, for times larger than Tw (N), the correlation (|43|) reaches 
a plateau. On larger time scales, the correlation (|^) should decay to its equilibrium value with a 

(2) 

cross-over time t w (N), corresponding for instance to the transition from diamond to graphite. The 

(2) 

time scale Tw (N) should increase with the number of sites, so that, with a macroscopic number of 
sites, no transition from diamond to graphite is observable on experimental time scales. On the other 
hand, t£\n) may remain small even in the thermodynamic limit, signaling the absence of ageing in 

(2) 

the graphite phase. In such a situation where the cross-over time scale t w (N) is not observable, the 
only relevant time scale is Tw(N). In order to calculate Tw (N) one should not use the prescription 
(44) but rather substract to the correlation ( f43| ) its plateau value for Tw\n) <C t <C t%(N). 



4 Brownian motion 

Brownian motion was shown in |l^] to be an ageing phenomenon. This model is simple and thus 
provides a simple test of the criterium of section |3|. Brownian motion in an euclidian space of 
dimension d has the following Glauber dynamics: Tj j = — 1/tq and T{j = l/(2dro) if i and j are 
nearest neighbors on a hypercubic lattice, tq is the microscopic hopping time. In this case, the 
Glauber matrix is symmetric and, because of translational invariance, its eigenstates are plane waves 

l*k> = -/=£ e * k - x l**>, ( 47 ) 
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where k = (nj, nd)2ir/L (we have assumed cyclic boundary conditions with a lattice spacing a). 
The matrix elements of the Glauber matrix on the plane wave basis are 



(48) 



where a runs over the d directions and k a is the component of the wave vector k on the a direction. 
The long time behavior of the random walk is readily obtained by expanding the cosine up to the 
second order, and the probability to find the brownian particle at time t w on site x is, in the long 
time limit: 

where the particle is at the origin at time t = (out of equilibrium state). As noticed first in p2| , 
the two times correlation functions of the position are proportional to the waiting time: 

t a 2 

CW(t w + t,t w ) = (x(t w + t).x(t w )) = (x 2 (^)) = (50) 

dT 

in the large t w limit. In order to make contact with [12|, one should adopt the scaling tq ~ a 2 in the 
thermodynamic/long time limits, which allows to recover the Langevin dynamics results and thus 
the existence of ageing in any dimension for the position-position correlation function. 

In order to make contact between the correlation (^) and the general form of the correlation 



(f43[), we should choose the overlap g(x) = x . With this choice of the overlap, (43) and (pC|) are 
both related to the same physical quantity: the mean square displacement. Using this specific form 
for the overlap, one can calculate Coo(t), integrate with respect to t and show that 7^(00) = 00, in 
agreement with the existence of ageing in this system for the position correlations. 
One could also be interested in the following correlation: 

C^(t w + t,t w ) = (5^ tw+t)fi 5 x{twh0 ) =V(0,t w \0,0)V{0,t w + t\0,t w ). (51) 

In the large t, t w limit, this correlation is found to be 

cP > (( „ +MroH (^_J_^. m 

In the thernodynamic limit, a — > 0, tq ~ a 2 , so that is vanishing in the thermodynamic limit, 
signaling the absence of ageing for the correlation (|51~p. The overlap associated to the correlation 
( |5TD is the local one: q x = 5 Xj o- We can thus calculate t w (oo) using the integration d46| ) over the 
density of eigenvalues. In order to calculate the behavior of t w (N), we consider a finite system of 
volume V = L d with a lattice spacing a. The number iV of sites is iV = (L/a) d . By thermodynamic 
limit, we mean that L is fixed, a — > and tq ~ a 2 . Using (|46|), we have then 
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where the summation is restricted over the wave vectors k such as 1/L < \k a \ < n/a. Clearly, the 
dominant contribution to (|53|) comes from the smallest eigenvalues only. It is thus legitimate to 
expand the eigenvalues up to the quadratic order in k: 

(54) 



One then obtains 



If d 



^w~^b^ 2 C^' (55) 



rW(JV)~»r„(ln(£)-ln(i)). (56) 
Since to ~ a 2 , Tw (N) — > in the thermodynamic limit in two dimensions. Now if d ^ 2, 

Since to ~ a 2 , this expression goes to zero in the thermodynamic limit whatever the dimension. 

We have thus shown that, in the case of random walks, and in the case of the correlations (^0|) 
and (|5l|), the existence/absence of ageing in the thermodynamic limit is consistent with the general 
argument of section |3|. 

5 The BM model at zero temperature 
5.1 Diagonalization of the Glauber dynamics 

We assume that the traps have been ordered such as E\ < E2 < ... < En. Unlike the REM-like trap 
model, we do not make any asumption about the trap distribution. The finite temperature dynamics 
is 

Ti ' j = N 1 + exp -£,)) (58) 
if i 7^ j and satisfies detailed balance. The zero temperature form of the Glauber matrix of this 
model is very simple since Z^j = — i + 1, Tj j = if i > j and Tij = 1 if i < j at zero temperature. 
The spectrum is {0,-1/ N, ...,(— N + l)/iV}. The eigenvector |*0o) for the eigenvalue A = is |ei). 
It is easy to show that the eigenvector \ifi-k) associated to the eigenvalue A = —k/N is 

l^-fc}=-^|efc} + -y=|e fc+ i). (59) 

We can now express the pure states in terms of eigenvectors 

|ei) = |Vo) (60) 
\e k ) = |^o) + V2 (|^-i> + ... + • (61) 
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This means that, at zero temperature, a pure state \et) [k is the /c-ieth traps, with E\ < ... < En] 
overlaps only with eigenstates associated to eigenvalues 0, — 1/N, (— k + 1)/N, and this overlap 
does not depend on the eigenstate. Compared to the result of section 2.2.4 the structure of the 



decomposition of a pure state into eigenstates of the Glauber dynamics is very different from the 
REM-like trap model. 

If one starts from a pure state | (0)) at time t = 0, the state of the system at time t is given by 



fc-i 



\e k (t)) = (1 - e-*/")|e x ) + £ (eH+W* _ e -lt/N^ ^ + e (-frH)*/tf | ejfc) . (62 ) 

1=2 

We recover the fact that this model can reach its ground state in a time of order N at zero temper- 
ature, whereas in the limit of zero temperatures, and for the REM-like trap model, the time spent 
in a trap i is proportional to exp — /3-Ei, which goes to infinity if /3 — ^ and N fixed. 

5.2 Two times autocorrelation functions 

The two times autocorrelation functions were calculated in [O] by means of a Laplace transform 
formalism. The result is 

HfotuO-r^- (63) 



if 1 <C t, t w . We rederive this result in the Appendix using a treatment different from [15], and based 



only on the diagonalization of the Glauber matrix. The BM dynamics thus exhibits ageing. 
5.3 Relation to the general argument 

Provided only self overlaps are non zero (qij = Sij), the correlation function introduced in section 
H is Coo(t) = j, provided t <C N. The typical time t^ m (oo) below which ageing occurs is thus 
logarithmically divergent. On the other hand, this can be also shown directly from the spectrum: 

JV-l 

.BM 



(AO = £ i (64) 



i=i 

which diverges logarithmically. 

6 One dimensional Ising chain 

As far as the finite temperature behavior of the two times autocorrelation function is concerned, this 
model exhibits interrupted ageing (?J since at any finite temperature the correlation length is finite, 
of the order of £ = exp (2J/T), ageing does not exist after a time independent of the system size. It 
is thus interesting to apply to this case the general argument of section ||. We first begin with the 
definition of the Glauber dynamics. 
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6.1 Glauber dynamics of the one— dimensional Ising chain 

We use the dynamics introduced in The master equation of this single spin flip dynamics is 
d ( N \ N 

^(Mi*) = - I X^M) j KM)*) + -,^N})p({cri, -,-o-i, ..,a N },t). (65) 

\i=l / i=l 

As a consequence of detailed balance, 

m(W}) = - (1 - (Ji tanh (/? J(a i+ i + <Tt-l)) • (66) 
Notice that, unlike the case of the REM-like trap or BM models, there is no need to multiply the 



matrix elements by 1/N. It is shown in [21] that in the infinite temperature limit, the eigenvalues of 
the Glauber matrix are negative integers between — N and 0, with degeneracies given by the binomial 
coefficients. The largest relaxation time in the infinite temperature limit is equal to unity whatever 
the system size, and there is thus no need to rescale the coefficients of the Glauber matrix. 

6.2 Relation to the general argument 

In the case of a spin model, the overlap q({cr}, {&'}) is given by 

1 N 

Q(WhM) = W J2w'i- ( 67 ) 

i=l 

Since ageing is interrupted, we should observe a saturation of the time T^f m9lD (N) as N increases. 
Moreover, Tl sin9lD {N) should increase as the temperature decreases. In practise, we work with 
an open chain and we take advantage of the existence of the global Z2 symmetry (invariance of the 
Glauber dynamics under the {a} — > {— a} transformation in a zero magnetic field), and the reflection 
symmetry of the chain. The Glauber matrix is diagonalized in each of these four symmetry sectors. 
The results are plotted on figure |. The difference rl sin9lD {N) - Tl sin9lD (N - 1) decreases as N 
increases, in agreement with the existence of interrupted ageing. 

7 SK model 

This model was proposed in 1975 as an "exactly solvable" model of spin glasses For reviews on 
this model and its developments, we refer the reader to the reviews f23j. The Hamilt onian of the SK 
model is 

// Y^.l, J a,a J . (68) 

where the sum is carried out over all the pairs of sites (infinite dimensional model) and the quenched 
random interactions are distributed according to a gaussian distribution 

2 \ 



P(J M .)=(^) 1/2 exp(-^j. (69) 
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This model is a spin glass if T < 1 and a paramagnet otherwise. The two times correlation function 
were shown in || to exhibit ageing in the glassy phase. We are going to apply our general argument 
of section ||| to the case of the REM-like trap model. To do so, we calculate numerically t^ k (N) for 
a large number of disorder realizations. The result is plotted on figure |(| as well as the differences 
A at = t w (N + 1) — t w (N) for N = 3, 4, 5, 6 sites. We observe that A^v increases at low temperatures 
and decreases at high temperature. The crossing point on the insert of figure |6| should be identified 
to the glass temperature. This temperature is not equal to unity since the system sizes are very 
small. Notice that even though the sizes are quite small, a large number of disorder configurations is 
required to have a significant statistics (of the order of 500000 independent disorder configurations). 
The computer time required was equivalent to 1000 hours of CPU time on a single processor. Our 
results are thus consistent with the existence of ageing in the low temperature phase, and no ageing 
in the high temperature phase. 



8 Conclusion 

We have thus carried out a detailed analysis of ageing from a Glauber dynamics point of view, with 
an emphasis on finite size effects. We have shown how to extract information from the Glauber 
matrix about ageing, and how to analyze finite size effects. In some simple cases with only self 
overlaps, the behavior of the time after which there is no ageing only involves the knowledge of the 
average density of eigenvalues pjv(A) when A — > - . This is the case for the REM-like trap model, 
in which case we have carried out a careful analysis of the density of small eigenvalues. Our analysis 
was confirmed by numerical diagonalizations of Glauber matrices of small clusters. In this case, we 
were able to related the typical behavior of the density of eigenvalues to the spectral dimension of 
the corresponding diffusion process. 
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Franz and N. Wingreen. The authors acknowledge B. Doucot for his interest to this work and his 
encouragements. R.M. also acknowledges the hospitality of NEC Research Institute at Princeton 
where part of this work was done. 



9 Appendix 

We are going to calculate the two times autocorrelation function II(t, t w ) using a method based only 
on the exact diagonalization of the Glauber matrix at a zero temperature. Starting from an infinite 
temperature state at time t = 0, 

1 N 

i*(o)> = ^Ei e *>. ( 7 °) 
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with \a) a pure state, the natural definition of the two times autocorrelation function within this 
model is 

N 

U(t,t w ) = J2'P(i,t w )expA i t, (71) 

8=1 

where V(i,t w ) is the probability to find the system in the trap number i at time t w and Aj is the 
escape rate from trap i, that is the diagonal element of the Glauber matrix Aj = X^. This definition 
of the waiting time dependent autocorrelation function is similar to the one used in Using the 
expression (|59|) of the eigenvectors and the decomposition (|60|) and (|6l] ) of the pure states into 
eigenvectors, we can calculate in a straightforward fashion the evolution of the system, starting from 
the state fl7|): 

\*(t w )) = 1 - (1 - ^)e- l -' N |ei) + - ]T e~ kt -l N \{N - k + l)e^ N - (N - k)] \e k ). (72) 

^ ' k=2 

The probabilities V(i,t w ) are easily recovered: V(i,t w ) = (ei\^(t w )). After straightforward calcula- 
tions, the dominant behavior of U(t,t w ) is obtained, in the limit 1 <C t,t w <^ N: 

W,t w )~-^-, (73) 

which is equivalent to the result of WjQ- Notice that if t w is of the order of the number of traps, the 
system has reached equilibrium: the dynamics is stationary (the system is stuck in the lowest trap) 
and there is no ageing. 
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Figure captions 



Figure 1: 

Density of eigenvalues in the Arrhenius regime for the REM-like trap model and fit to the expo- 
nential form (/3 = 9,iV = 75 traps). The insert shows the same quantity, but with also the smallest 
eigenvalues. For these small eigenvalues, the Arrhenius picture is not valid. 



Figure 2: 

Density of small eigenvalues of the models of section |2|, for (3 = 5. (a): infinite dimensional case 
(REM-like trap model); (b): one dimensional case, with "disordered" traps; (c): one dimensional 
case with "ordered" traps. The lines in cases (b) and (c) are parallel, and in good agreement with 
the Arrhenius law, whereas we observe a significant deviation to the Arrhenius law in (b) case. 



Figure 3: 

Exponent a deduced from the fit (^) of Pb (In |A|). The circles correspond to the infinite dimensional 
case (REM-like trap model), the squares to the one dimensional "disordered" case, and the diamonds 
to the "ordered" one dimensional case. The solid line corresponds the the Arrhenius regime (in other 
words, to site-disordered random walks above the upper critical dimension 2), and the dashed line 
corresponds to the one dimensional renormalized Arrhenius regime (in other words, to site-disordered 
random walks in one dimension). 



Figure 4: 

Structure of the eigenstates of the REM-like trap model for 20 traps. The label i is an eigenvalue 
label, with eigenvalues in growing order and j labels the traps (the depth increases if j increases) . 
(a) represents for the 20 pure states, (b) represents the depth of the traps, (c) represents 

the 20 eigenvalues between —1 and for 4 different temperatures (circles: (3 = 2, squares: (3 = 5, 
diamonds: (3 = 10, triangles: (3 = 20). 



Figure 5: 

Variations of rlf in9lD (N) as a function of the inverse temperature (3 for N = 3,4,5,6,7,8,9 sites. 
The insert shows the differences rl sin9lD (N) - Tl sin 3 1D (N - 1) for N = 4,5,6,7,8,9 sites (curves 
a,b,c,d,e,f respectively). 



Figure 6: 

Variations of t^ k (N) as a function of the inverse temperature (3 for a N = 3, 4, 5, 6 SK cluster. The 
insert shows t^ k (N) - t^ k (N - 1) [a: iV = 6, b: N = 5 c: N = 4]. 
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